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Imidacloprid, one of the most commonly used insecticides, is highly toxic to bees and other beneficial 
insects. The regulatory challenge to determine safe levels of residual pesticides can benefit from information 
about the time-dependent toxicity of this chemical. Using published toxicity data for imidacloprid for 
several insect species, we construct time-to-lethal-effect toxicity plots and fit temporal power-law scaling 
curves to the data. The level of toxic exposure that results in 50% mortality after time t is found to scale as t 1 7 
for ants, from t 1 6 to t 5 for honeybees, and from t 1 46 to t 2 9 for termites. We present a simple toxicological 
model that can explain t 2 scaling. Extrapolating the toxicity scaling for honeybees to the lifespan of winter 
bees suggests that imidacloprid in honey at 0.25 fig/kg would be lethal to a large proportion of bees nearing 
the end of their life. 



The European Commission recently voted to place a moratorium on the use of three neonicotinoid insecti- 
cides (i.e. imidacloprid, thiamethoxam and clothianidin) after the European Food Safety Authority (EFSA) 
determined that there were "high acute risks" to bees from these products through several exposure routes 1 , 
e.g. by direct dust exposure during seed planting operations 2 , or by ingestion of residues in guttation droplets 3 , 
nectar and pollen from treated plants 4 . In their assessment of the risks from imidacloprid on bees, the EFSA 
considered toxic endpoints for acute (3.7 ng/bee) and chronic exposure (20 ug/kg or L, henceforth ppb) that 
represent a rough consensus of the toxicity studies reported in the literature (see review in 5 ); however, they 
cautioned that there are no guidelines for chronic and sublethal exposure testing in bees, and expressed concern 
regarding the uncertainty about the biological significance of such exposures. 

One way of resolving that uncertainty is to consider time as one of the exposure factors, such as concentration. 
Indeed, toxicity scaling with time can lend insight into the effects of residual pesticide contamination and can 
provide guidance to regulators for determining acceptable contamination levels. 

Neonicotinoid insecticides are based on the natural toxin nicotine 6 , and are of particular concern because they 
bind virtually irreversibly to the nicotinic- acetylcholine receptors in the insect's nervous system 7 , so the damage 
can accumulate, and therefore the toxic effects can be reinforced with chronic exposure 8 . Time-to-effect studies 
lend themselves to a simple time- dependent power-law empirical model which can guide expectations for field 
toxicity effects 9 . Here we attempt to quantify risks to honeybees by considering the time- dependent toxicity of 
imidacloprid, the most ubiquitous of the neonicotinoid pesticides and the number one insecticide sales in the 
world 10 . The lethal time-to-effect relationship can be expressed as a power law 



LT50xD t p 



(1) 



where the toxic effect is the time when half of the insects succumb to the toxin (LT50), D is the dose of toxin per 
unit time, and P is the power law exponent. A scaling relationship, such as Equation 1, is useful to extend predicted 
effects to low doses where experimental data do not exist. 

There are only a few studies with honeybees that follow the insects beyond ten days, e.g. 11 , so we also looked at 
studies with ants and termites to corroborate the temporal toxicity dependence that we find in the honeybee 
studies. Finally, we consider a simple toxicological model to explain the experimental data. 
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Time-to-effect studies with terrestrial insects 

The studies analyzed here include cases where the exposure was by 
direct ingestion of the pesticidal toxin by adult honeybees and ants, a 
study where the toxin was fed to larvae of stingless bees, and finally a 
study of toxic exposure through treated soil for termites. Rust et al. 12 
reported time- dependent toxicity measurements for imidacloprid on 
Argentine ants, Linepithema humile (Mayr). Data for termites 
(Reticulitermes flavipes (Kollar)) were taken from Ramakrishnan 
et al. 13 , who reported mortality and feeding inhibition for termites 
observed up to 21 days in different soils treated with imidacloprid. 
The stingless bee Melipona quadrifasciata anthidioides Lepeletier, 
provides a good subject for studying the effect of the toxin on devel- 
oping insects because the larva feeds only from food sealed inside the 
cell, so the total dose ingested can be accurately determined. The 
study by Tome et al. 14 followed these bees for 45 days from the egg 
until just after emerging as an adult over a wide range of toxic doses of 
imidacloprid. 

For honeybees {Apis melifera L.), a paper by Suchail et al. 15 and 
another by Dechaume-Moncharmont et al. 16 report time-to-effect 
mortality data after chronic exposure to imidacloprid. However, 
there has been controversy about Suchail et al. results because two 
other studies did not show such a high sensitivity to the toxin 1718 . 
Hence, we also include time-to-effect measurements from a study by 
the UK Department for Environment, Food and Rural Affairs 
(DEFRA 19 ), and a recent compilation of many researchers' acute 
median lethal doses (LD50) values for 24, 48, 72, and 96 h periods 
(FERA 20 ). 

Case studies 

Poisoning of adult insects by ingestion of imidacloprid. In order to 
plot the acute toxicity data along with chronic data on the same graph, 
average daily dose rates were calculated by dividing the acute LD50 
values by the time intervals that were used for the measurement. For 
the chronic data, we use the reported concentration and consump- 
tion rates to obtain the daily dose rate. In both cases, the reported 
consumption of toxin (ng/day) is plotted versus the time of exposure 
until half the bees have died (LT50). 

Most studies follow the insects for 10 days or less, so at low toxin 
concentrations the LT50 time cannot be reached. This is true for 
several of the data sets we analysed 20 , so whenever mortality is 
reached between 15% and 50% with respect to the control, we extra- 
polate from the last two time points to estimate the LT50. The paper 
by Dechaume-Moncharmont et al. 16 contains a couple of data points 
with time-to-effect numbers for relatively low sublethal concentra- 
tions. Data from their paper (Figure 1) show that 4-8 ug/L (ppb) 
concentrations of imidacloprid in syrup eventually kill 50% of bees in 
about 30 days. Figure 1 also shows how we determine the LT50 point 
from the mortality time series using Abbott's formula 21 . LT50 times 
falling between reported sample points are linearly interpolated from 
neighboring points. 

The data for poisoning of adult insects after ingesting imidacloprid 
are plotted in Figure 2. A power law best-fit line for each data set is 
also shown. The slope of the fitting line is — 1/P, the exponent of 
power law Equation 1 . For example, the LT50 times for the Argentine 
ants fit remarkably well a simple power law model with roughly t L? 
dependence (r 2 = 0.97, p = 0.013). This adds confidence to the utility 
of the empirical power law model for imidacloprid with these social 
insects. Recall that simple accumulation to a toxic threshold (Haber's 
rule) would appear as directly proportional to time (t 10 ) 22 . An expo- 
nent larger than one can be interpreted as coming from damaging 
secondary physiological effects that develop over time. We will come 
back to this interpretation later. 

The DEFRA honeybee data stand out as not falling within the 
averages of many other researchers in the FERA report 20 ; about a 
tenfold difference, despite a good fit to the power model (r 2 = 0.99, p 
< 0.001). There appears to be a wide variation between experiments, 
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Figure 1 | Time-to-effect curves for imidacloprid fed to honeybees at 4 
and 8 ug/L in syrup (adapted from Dechaume-Moncharmont et al. 16 ). 

which could be due to the individual colony or strain of bees' sens- 
itivity to this insecticide 23 ,or by other factors not considered in the 
experimental design (temperature, season, condition of bees, nutri- 
tive status, etc.) The DEFRA study reports time-to-effect numbers 
and consumption data, so it provides a good data set for determining 
time scaling even if the DEFRA bees are not particularly represent- 
ative of most other bees tested. The mortality for the DEFRA bees can 
be very well explained by the power law formulation with about t L6 
dependence, very similar time exponent to the ants. Finally, the 
average LD50 values for 24, 48, 72, and 96 hours, compiled from 
many researchers, also fits the power law model (r 2 = 0.98, p < 
0.001), here about f. The Dechaume-Moncharmont et al. 16 data 
are included in this power fit line. 

Dechaume-Moncharmont et al. 16 results, using concentrations of 
4 and 8 ppb, (0.08 and 0.16 ng/day, respectively) are relevant to 
field- realistic exposure levels. Residue concentrations of imidaclo- 
prid in pollen and nectar from treated agricultural crops are fre- 
quently reported in the 1 to 30 ppb range in France 24 and the US 25 . 
Bees foraging on those crops could be subject to the test levels used by 
Dechaume-Moncharmont et al. 16 , which caused 50% mortality in 
about 30 days. 

The Suchail et al. 15 results show a much higher sensitivity to imi- 
dacloprid than the trends we see with the other honeybee data sets in 
Figure 2. Rather than a t 2 dependence, lethal effects seem to scale as 
toxic concentration times f (r 2 = 0.87, p = 0.021). It is tempting to 
disregard the Suchail et al. 15 results; however one test site in trials 
reported by Schmuck 18 also showed very high sensitivity to a meta- 
bolite of imidacloprid (6-chloro-nicotinic acid) as well. That study 
reports 10-day mortality at the Germany II location: control 10%; 
0.004 ng/bee/d 67%; 0.038 ng/bee/day 77%; 0.38 ng/bee/day 97%. 
Other locations did not show such high sensitivity to the insecticide. 
It is worth considering secondary stressors that could result in a 
higher sensitivity to the toxin. Suspects would be bacterial or viral 
pathogens or the Nosema microsporidian 26 , or even synergistic inter- 
action with other insecticides 27 . 

Compelling new research by Di Prisco, et al. found that very low 
doses of both the neonicotinoids clothianidin and imidacloprid 
caused immune suppression in honeybees as evidenced by replica- 
tion of the deformed wing virus (DWV) 28 . They provided evidence 
that the neonicotinoids were downgrading the innate immunity 
pathway governed by NF-kB. Replication of the virus was found to 
increase with both dose and time. 

Pathogen interactions with a host could lead to several time- 
dependent processes. First, damage to the host immune system has 
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Figure 2 | Time-dependent toxicity from ingestion of imidacloprid for honeybees and Argentine ants. Published insect mortality results from several 
researchers plotted as time-to-effect versus daily dose per bee. Field exposure range assumes <5 ppb a.i. at 20 ul consumption/ day. 



to occur for a pathogen to get a foothold. Then, before the host 
succumbs to the pathogen, the infectious organism must grow and 
multiply to lethal levels. Damage from the pathogen itself may take 
time to manifest in the host organism. Several of these time- 
dependent processes, initiated by pesticide exposure, could produce 
a high order time dependency like that observed by Suchail et al 15 . 

Addressing this question, the interaction between Nosema and 
imidacloprid was studied by Alaux et al. 29 . The study fed 200,000 
Nosema spores to bees to initiate infection and then subjected the 
bees to various doses of imidacloprid. There was some interaction 
between the pesticide and the pathogen, but the study only fol- 
lowed the bees for ten days, so it is hard to draw conclusions at 
field- realistic imidacloprid doses, as the Nosema infection domi- 
nated the experiment at this short time scale. Nevertheless, two 
data points from that study are plotted in Figure 2 (crosses) where 
estimates could be made of the LT50 time with the addition of 
imidacloprid, using the Nosema infected bees as the control. 
However, imidacloprid- alone concentrations from 0.7 to 70 ppb 
in that study resulted at day 10 in only about 7% to 12% mortality 
above the control, not sufficient to render an LT50 estimate. 

Poisoning by ingestion of toxin in larvae. When imidacloprid is fed 
to developing larvae, the dose-time dependency is markedly 
different. Stingless Melipona bee workers will succumb to the 
toxin's effects between day 11 and day 19 of their development for 
a toxic dose range spanning more than two orders of magnitude: 
from 14 to 2800 ng/bee, as shown in Figure 3 14 . 

After the larval stage, the survival probability for the individuals 
making it to the pupae stage improves dramatically. At the lowest 
dose tested, 5.6 ng imidacloprid a.i. in 130 |il of larva food before 
sealing the cell, more than half of the bees survive well into adult- 
hood, yet a third of that cohort died between day 1 1 and pupation at 
day 23. In this case the authors observed an almost dose-independent 
effect that caused mortality during a specific period of larval develop- 
ment. The time dependence for this section of the toxicity curve is 
t 13 \ (r 2 = 0.80, p < 0.001). This might also be (at least partially) due to 
the changing physiology of the developing insect, rather than to the 
cumulative effect. Moreover, the food uptake increases in the grow- 
ing larva, thus increasing the exposure. 



It is worth noting that had the authors followed the mortality of 
the larva for just ten days, as is common with adult honeybee toxicity 
studies, they would have missed the main feature of delayed toxicity 
which started at day 11. Similarly, mortality of honeybees in the 
Dechaume-Moncharmont et al. 16 study started to occur about day 
19 of exposure in the 4 ppb concentration (Figure 1). 

The greater survivability once the bees reach the pupae stage can 
be explained by two factors. First, the pupae are no longer consuming 
the toxin-laced larval food supply, so the internal concentration of 
toxicant would plateau or even decrease. More importantly, the 
insects continue to develop neurologically. A growth spurt of syn- 
aptic dendrites is known to occur at metamorphic transitions 30 which 
would allow repair and recovery from previous damage in the sur- 
viving cohort. This mechanism may mean that larval sublethal expo- 
sure is only relevant for that stage of development, but largely 
irrelevant for insects that are further exposed as adults. 



100 
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Figure 3 | Time-dependent toxicity for larval stingless bees fed 
imidacloprid (adapted from Tome et al. 14 ). 
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Poisoning of termites in imidacloprid-treated soils. The previous 
studies refer to oral ingestion of imidacloprid. Moreover, time-to- 
death data for Reticulitermes termites exposed by contact to several 
imidacloprid-treated soils reveal a similar toxicity pattern (Figure 4). 
In this case, the power law model renders exponent values that range 
from 1.46 to 1.62 for loam and sandy-loam soils, respectively, to 2.6 
in sand soil and 2.9 in clay-loam soils (r 2 range 0.73 to 0.94, p range 
0.005 to 0.045). The authors of that study 13 explain that the different 
adsorption properties of each soil tested can account for the different 
exposure of the termites, hence the variable patterns of chronic 
toxicity. Besides the lethal endpoint, they also report significant feed- 
ing inhibition for termites in the imidacloprid-laced soils, possibly a 
contributor to the lethal effect, as observed also with bumble bees 31 . 
Although wide-ranging in absolute toxicity depending on the soil 
type, the time dependence of the toxic effects is similar to the other 
adult insects we have looked at above, with time exponents near 
two. 

A Toxicity Model for neonicotinoids. Imidacloprid and all other 
neonicotinoids are designed to bind to the nicotinic acetylcholine 
receptors (nAChRs), which are located on the dendritic terminal 
membrane at the synapses of the insect nervous system 6 . When the 
neonicotinoid molecule binds to the nAChR the ion channel opens, 
as it would if it was stimulated normally by the acetylcholine (ACh) 
neurotransmitter; in other words, neonicotinoids are agonists that 
compete with the natural neurotransmitters. The neonicotinoid 
molecule remains bound to the nAChR in insects, holding the 
channel open, whereas normally ACh remains bound for only 
about a millisecond 32 before dissociating from the receptor and 
being removed from the synaptic junction by the enzyme 
acetylcholinesterase (AChE). In mammals and other vertebrates, 
the lesser affinity of neonicotinoids for their nAChRs appears to be 
related to the different configuration of the subunits that make up 
this receptor, so the toxin binding is weak and/or does not last as long 
as in insects 6 . 

A toxicokinetic-toxicodynamic model (TKTD) to explain the 
time-cumulative toxicity of imidacloprid shown above is proposed 
here. The model is also applicable for all neonicotinoids. The tox- 
icokinetic portion of our model consists of a rate equation that 
describes accumulation of bound toxin on the synapses. We take 
the toxic exposure dose, D(t), to be either a single initial dose or a 



continuous, much lower, chronic dose. The change in total-body 
toxic load, C, can be written 



dC 
~dt' 



--D(t)-C/T M + C B /T D 



(2) 



where x M is the metabolic decay time. Neurotoxins produce an effect 
only when bound to specific receptors, so the toxicodynamic part of 
the model relies on the relative velocities of binding and dissociation 
to and from the receptor. Therefore, the amount of bound toxin, C B , 
will depend upon the total toxic load and the time constants for 
receptor binding, t a , and dissociation, t d . 



dC B 
dt 



--C/t a — C b /t d 



(3) 



We assume here that there are plenty of binding sites and that the 
binding rate does not decrease because a large fraction of the sites 
become filled. However, one would expect that a very small percent- 
age of the ion channels bound with toxin and held open continuously 
would be damaging to the organism; by way of analogy, it is like 
having a power switch turned on for a long time, squandering energy 
unnecessarily. Desensitisation, endocytosis or neuronal adaptation 
to the bound neonicotinoid could also occur, which we do not 
account for explicitly in our model other than through the empir- 
ically determined time constants. 

It is less clear how the amount of bound neurotoxin produces 
toxic effect over time, even if there is ample evidence of it occur- 
ring 33 . An explanation of the mechanisms leading from sublethal 
effects to delayed lethal effects over time needs to consider the 
specific mode of action of the toxin 22 . A likely candidate is the 
phenomenon of excitotoxicity, where prolonged and excessive 
stimulation of nerve cells leads to cell death from a variety of 
mechanisms 34 . For example, continuous muscle stimulation by 
the affected nerve can lead to death of the muscle cell. Excessive 
stimulation of the synapse can lead to Ca ++ imbalance and demise 
of the nerve cell itself. Both of these mechanisms are the result of 
the accumulation of insults to the affected cells from synaptic 
stimulation over time. We can express the damaging biological 
effect, E(t), as 



E(t) oc C B dt 



(4) 
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Figure 4 | Toxicity of imidacloprid-treated soil to termites (data from Ramakrishnan et al. 13 ). 
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Figure 5 | Model results for chronic exposure of honeybees to 
imidacloprid at 0.16 ng/day. Unbound total-body toxic load (solid line) is 
a small fraction of the daily intake because most of the toxin is metabolized 
while the amount of bound toxin (short dashes) grows linearly with time. 
Relative biological damage grows as t 2 . The model is scaled such that 
biological damage is unity at the point the LT50 is reached for honeybees. 

We interpret a measured toxic endpoint to be where E(t) reaches 
an arbitrary threshold, e.g. E(t) = 1.0 corresponds to LT50. 

To apply the model, we need estimates of kinetic time constants 
?d> t m and t a . Estimates for x M and x A are obtained from the radio- 
carbon metabolic analysis done by Suchail et al. 35 . Although the 
parent compound is quickly metabolized and excreted with a lifetime 
of about 5 hours, several toxic metabolites remain with a much 
slower apparent decay time 35 . Approximately 10% of the initial dose 
is still found in the insects after 72 hours as metabolites. This implies 
that the ratio t m /t a ~ 0.1. An accurate estimate of x D was not 
possible, because the bees were only followed for 72 hours in 
Suchail et al. 33 experiments, not long enough to establish the decay 
rate of the bound component. For irreversible binding, x D would 
approach infinity. 

An example of this simple model is shown in Figure 5 for chronic 
exposure of honeybees to imidacloprid at environmentally relevant 
oral doses of 0.16 ng/day. The model is scaled so that at the observed 
LT50 of honeybees, the model's biological damaging effect, E(t), is 
approximately equal to one. In the chronic case, the total body toxic 
load remains at about 20% of the daily dose because the toxin is being 
continuously ingested and rapidly metabolized. However, the bound 
toxin that is trapped on synapses grows linearly, while the biological 
damage increases with t 2 , and this outcome of the model matches 
very well the mortality observed in terrestrial insects (see Figures 2 
and 4). 

Further considerations of the model utility are discussed in the 
Supplementary Information (SI). 

Discussion 

At present, regulatory guidelines primarily assess the survival of 
adult honeybees after a short exposure to pesticides, acutely up to 
4 days, and with chronic caged bee studies up to 10 days 36 . However, 
some authors have stressed the importance of also testing for longer 
chronic toxicity, larval toxicity, and sublethal effects of pesticides 37 . 
Without explicit long term tests, we are left to estimate long term 
effects by extending the toxicity scaling that we find in shorter dura- 
tion studies. We can extrapolate the curve in Figure 2, beyond the 
Dechaume-Moncharmont et al. 16 data points with LT50 at about 30 
days, to the lifespan of wintering bees at 150 days 38 . The extrapolation 
suggests that daily ingestion of about 0.005 ng/day of imidacloprid 
would produce LT50 in 150 days. For bees consuming 0.02 g honey 
per day 4 , this implies a concentration of 0.25 ppb in honey as the 



lowest concentration capable of causing long-term mortality. Such a 
low concentration is about the detectable limit of this insecticide by 
LC/MS and below its limit of quantitation 24 . Hence, even with 
healthy bees, exposure to modest field- realistic residues of imidaclo- 
prid in pollen (range 0.5-30 ppb) and honey (range 0.7-13 ppb) 
could easily cause problems for summer bees and especially for 
longer-lived bees going through the winter. Certainly, not all pollen 
or honey are contaminated with imidacloprid, so the oral exposure of 
bees would be scaled down by a certain factor, since the average 
frequency of detection in pollen is between 9% and 40% for the US 
and France, respectively 24 ' 25,39 . 

Extrapolation of the toxicity curve should not be considered defin- 
itive, but rather cautionary. The toxicity scaling observed up to 30 
days may not continue, as there is uncertainty in the actual decay rate 
of bound toxin molecules. Our model results suggest that decay times 
as short as 10 days are consistent with the observed data (See Figure 
S2). It is also likely that there are differences in the rate of food 
consumption and physiological differences between summer and 
winter bees that could affect the mortality rate 40,41 . 

Colony health may not suffer significantly due to mortality of 
older bees as long as the colony remains otherwise healthy with a 
productive queen that is able to maintain the colony population. It is 
well known that colonies can compensate for the losses of worker 
bees. Indeed, full colony field studies on neonicotinoid treated crops 
show not much difference in performance between the insecticide 
treated and non-treated colonies 41 ' 42 . Also, a field study investigating 
imidacloprid-treated maize found the level of contamination in 
stored honey was between 0.05 ppb and 0.5 ppb, and yet mortality 
rates in apiaries were inversely correlated with the surface of maize 
fields treated but not with imidacloprid 43 . Additionally, risk analysis 
studies and reviews of neonicotinoid use have downplayed the 
importance of the neonicotinoids on honeybee losses, especially as 
the sole agent 517,44 . However, the explicit connection between innate 
immunity loss and the neonicotinoids 28 , was not known or consid- 
ered in these studies. The potential for accumulation of neonicoti- 
noids at receptor sites to the level where loss of innate immunity has 
been demonstrated may be much more likely for sublethal exposure 
to these chemicals than direct poisoning. Furthermore, it must be 
recognized that worker losses do occur, and declining homing rates 
can be accurately measured using radio frequency identification 
(RFID) tags 45 " 47 . 

From the time-to-effect studies we have examined, in all cases it 
appears that imidacloprid lethal effects accumulate in the insects. 
Observations by many authors show that the longer the exposure 
time, the less amount of total chemical needed to kill the insects; in 
other words, the LD50s decrease with exposure time, as shown in 
Figures 2 to 4. Whether it is enhanced, cumulative, or delayed tox- 
icity, all these terms describe this situation when the power law 
scaling exponent is greater than one. One explanation is that the 
damage at cellular/tissue level might occur immediately after the first 
administration of the toxin, even if no permanent binding is 
observed. However, this damage might produce no immediate visible 
effect - mortality of the insect will occur later on. Thus there might be 
acute poisoning, with no mortality, during the initial phase of the 
experiment and its delayed lethal effect would increase the exponent 
value even if no toxin accumulation occurs. 

Insecticides with delayed toxicity pose a particularly vexing prob- 
lem when attempting to protect beneficial insects. Separation of 
residual toxin in time and space between target insect pests and 
pollinators becomes increasingly important as the order of time 
dependence increases. For example, consider a case where applica- 
tion of the pesticide is designed to kill target insects within 2 days. 
Assume that the pollinator lifespan is 50 days and we wish to not 
harm the pollinator. Further, assume that the 48 hour LD50 is the 
same for the target insect and pollinator (since they are both insects). 
The usual approach to protecting pollinators is to apply insecticide 
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when bees are not present so the full force of the insecticide can act 
upon the target pest while only a residual amount will remain in the 
crop when bees encounter the pesticide. Table 1 shows the level of 
pesticide contamination compared to the application rate required to 
protect a long-lived pollinator, and the safety factors that need to be 
applied. 

Compounds that exhibit enhanced toxicity require much more 
care for the toxin to be well-eliminated from the pollinator's envir- 
onment. Systemic insecticides pose a particular risk in this regard, 
because once the toxin has been taken up in the plant, it can remain 
there for an extended period, including plant bloom, so their residues 
translocate to pollen and nectar 48 ' 49 . Achieving the protection ratio 
required for pollinator safety may be impossible for long-lived com- 
pounds. Current regulatory standards emphasize a toxic threshold 
that does not include any time dependence. Utilizing a simple toxic 
threshold to define safe residual levels based on acute LD50s, which 
works for many pesticides, has led to complacency when faced with 
compounds that exhibit enhanced toxicity, such as imidacloprid and 
related chemicals. Compounds with a t 2 toxicity scaling require pro- 
tection levels hundreds of time more stringent than those with 
merely a static lethal threshold (Table 1), because the eco- 
toxicological consequences of their delayed toxicity are potentially 
far-reaching. Many invertebrate species have long-lived stages of 
development in which delayed toxicity of residual pesticide could 
directly affect the organisms. 

We examined studies that addressed bee mortality after exposure 
to imidacloprid. Bees thus affected may not immediately die, but they 
may be practically useless to the colony, as they may be unable to 
forage or navigate 11,50 " 52 , and hence become lost and perish outside 
the hive. Naturally, while some bees are lost, others are produced to 
replace them, so the population in the colony is compensated, mak- 
ing it difficult to estimate the actual mortality rates. Quantifying 
behavioral effects is more difficult than counting dead bees in a hive, 
and may require the use of sophisticated equipment to track down 
each individual bee 53 . We would expect that time-to-effect scaling for 
behavioral effects would be similar to the mortality data, but would 
merely show up sooner. Extrapolating from individual bee mortality 
or neurological damage to whole-colony health remains an impor- 
tant unresolved question. Neurological damage can effect foraging 
success and hence colony nutrition 54 , or grooming behavior, learning 
and olfactory sensory abilities 50,55 , which are important colony- wide 
responses to parasites and pathogens. 

It is generally accepted that multiple pathogens ultimately bring 
down stressed colonies 56,57 . Cornman et al. 57 found that Colony 
Collapse Disorder (CCD) colonies were more likely to have higher 
levels of a wide variety of pathogens than weak, but non-CCD, col- 
onies. Not only were the levels of pathogens higher, but multiple 
agents were frequently found in combinations not typical of non- 
CCD colonies. It appears that the immune system in the CCD col- 
onies has gone awry, and some authors hypothesize this is the case 58 . 
Could it be a few parts per billion of insecticide that makes the 
difference? The Di Prisco study showed both imidacloprid and 
clothianidin adversely affect insect immunity and promote replica- 
tion of a viral pathogen, DWV, in honey bees at exposure levels 1 ppb 



or less with exposure of one to three days 28 . More research is still 
needed here, perhaps looking at pesticide interaction with KBV or 
AMPV virus, since these viruses show up in CCD colonies. 

Aufauvre et al. 59 looked at the interaction of the systemic insect- 
icide fipronil and Nosema infection, and found a strong synergy 
between the pesticide and the pathogen, similar to what happened 
with imidacloprid 29 . Bees with either Nosema alone or fipronil alone 
suffered mortality similar to the control, since low doses of toxin and 
Nosema spores were used. However, given enough time, the com- 
bination was especially lethal. The study followed the bees for 22 
days, the added time crucial for the delayed toxic and synergistic 
effects to show up. 

The choice of Nosema as the pathogen for pesticide-pathogen 
interaction studies is convenient because the infection is easily 
accomplished with a spore solution, and the progress of the infection 
can be followed with microscopic examination of the bees 26 . Viruses 
are much more difficult to use as the infectious agent because of the 
difficulty of diagnosing their presence and quantifying the infectious 
dose and infection progress. However, viruses are more ubiquitous in 
honeybee colonies than Nosema, and could easily have a similar 
synergistic interaction with pesticides on the bee's health. 

A final point to consider is the effect of imidacloprid on queen 
bees. Queen bees can live for several years. Queen bees also need to 
consume large amounts of food (e.g. honey and royal jelly) in order to 
churn out the quantity of eggs they lay every day. This makes queen 
bees especially susceptible to a cumulative toxin in the food supply. 
Queen failure was one of the precursors to colony mortality found by 
vanEngelsdorp et al. 60 . Equally, colonies of bumble bees (Bombus 
spp.) exposed to imidacloprid suffered an 85% reduction in produc- 
tion of new queens compared to control colonies 61 . The situation is 
aggravated when the bees are exposed to not just imidacloprid alone 
but two or more insecticides 27 , as it occurs in the real world. Queen 
failure, even if the bees succeed in raising a new one, will leave the 
colony without a fresh supply of young bees. A colony with pesticide 
stress may be relying on young bees to make up for those that dis- 
appear before their time. 

Conclusions 

Accurate time-to-effect scaling is an important tool for estimating 
the effect of chronic pesticide exposure to honeybees. In order to 
achieve this, chronic tests for pesticide toxicity to pollinators should 
be extended to 30 days or more and use time-to-effect measurements 
to clearly establish scaling. This is because in some experiments there 
were no deaths in either the control or exposed bees during the first 
10 days 16 , yet the toxic effect was clearly evident and different from 
the controls by day 30. Similarly, with stingless bee larva, the growing 
larvae required at least 1 1 days before the toxic effects of imidacloprid 
became apparent. The model shown here reproduces well the obser- 
vations made by many authors and explains the delayed and time- 
cumulative toxicity of imidacloprid based on its mode of action. 
Toxins that have enhanced time-dependent chronic toxicity, 
such as imidacloprid and other neonicotinoids, need to be well- 
characterized to avoid accidental poisoning of beneficial organisms. 
It follows that compounds with delayed toxicity should be avoided as 



Table 1 | Maximum residual insecticide contamination allowed to protect a long-lived pollinator (e.g. life-span 50 days) compared to the 
application rate when the insect pest is targeted in 2 days (Ratio pollinator/ pest = 50/2). Assumes identical sensitivity to the toxin for both 
target and non-target insects. Relative safety factors (x 3) required to protect the beneficial pollinator are indicated for each case 

Residual Concentration compared to 

Pesticide toxicity time treatment level that kills non-target 

dependence Description insects Include Safety Factor X 3 

f° Threshold level only- doesn't depend on time 1/1 1/3 

f 1 Accumulate to threshold with time - Haber's rule 1/25 1/75 

t 2 Enhanced or delayed toxicity 1/625 1/1875 
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pesticides because of the intrinsic difficulty they pose as envir- 
onmental contaminants. 
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